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^vq , Abstract We provide formulas for the moments of the real and complex non- 

central Wishart distributions of general degrees. The obtained formulas for the 
real and complex cases are described in terms of the undirected and directed 
graphs, respectively. By considering degenerate cases, we give explicit formu- 
^vq \ las for the moments of bivariate chi-square distributions and 2x2 Wishart 

distributions by enumerating the graphs. Noting that the Laguerre polynomi- 
als can be considered to be moments of a noncentral chi-square distributions 
^ | ■ formally, we demonstrate a combinatorial interpretation of the coefficients of 

\ the Laguerre polynomials. 

' Keywords Kibble's bivariate gamma distribution • Laguerre polynomial • 

Noncentral Stirling number of the first kind 



1 Introduction 

> 

[• — ' For t = 1, . . . , v, let Xt = (xu)i<i< p be a p-dimensional random column vector 

distributed independently according to the normal distribution N p (^ t , with 
mean vector /i t = (/iti)i<i<p and covariance matrix £ = (&ij)i<i,j<p- We de- 
fine the (real) noncentral Wishart distribution W p {v, S, A) by the distribution 



q 

of a p x p symmetric random matrix 



W = (wij)i<ij< p , Wjj = y^xuxtj, (1) 



t=l 
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where 

a = (5ij)i<i.j< P , — nuptj 

t=i 

is the mean square matrix. The distribution of W depends on /i t 's through A 
because its moment generating function is 

E[ e ^ ew )] = det(7 - 20E)S e *V-*eE)-*eA } (2) 



where is a p x p symmetric parameter matrix (jMuirhead I (|1982f> ). 

Note that A = if and only if fa = for all i. The Wishart distribu- 
tion with A = is referred to as the central Wishart distribution W p (u, S). 
Conventionally, the triplet (i/,E,f2) with J? = E~ 1 A is used for describing 
the noncentral Wishart distribution rather than (i>, S, A). The matrix i? is 
called the noncentrality matrix. In our paper, we adopt the triplet (y, S, A) 
for simplicity in describing theorems. 

For t = 1, . . . , z/, let ( v * J be a 2p-dimensional random column vector dis- 



Yt, 

tributed independently according to the normal distribution with mean vector 

I) and covariance matrix "/), where A and B arepxp symmetric and 

skew-symmetric matrices, respectively. The distribution of a complex-valued 
random vector Z t — (zti)i<i< p = X t + y— lYt is referred to as the complex 
normal distribution CN p (fi t , U) with mean /i f = {^u)i<i< P = £t + V - l^t 
and covariance matrix U = {uij)i<i.j< p — 1{A + \/ — lB). Actually, S is a 
"covariance" in the sense of 

<7ij = E[(z u - Vu)(ztj ~ fJ-tj)]- 

Here, the over line denotes the complex conjugate. From the complex random 
vectors Z t , we define the complex noncentral Wishart distribution CW p (is, S, A) 
as the distribution of a p x p Hermitian random matrix 

V 

W = (wij)i<ij< p , Wij = z tl zTj, (3) 

t=i 

where 

a = (5ij)i<i.j< P , — utiiHj 

t=i 

is the mean square parameter matrix. As in the real case, the distribution of 
W depends on /i t 's through A since its moment generating function is 

E[ e ^ ew )] = det(J - 0S)-"e t ^ I - es ^ l0A , (4) 



where Q is a p x p Hermitian parameter matrix. See 
central case. 



Goodman 
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The primary purpose of this paper is to obtain expressions for the moments 
E[w a bW c d ■ ■ ■ w e f] of the real and complex noncentral Wishart distributions in 
terms of graphs, where a,b,c,d,...,e,fG {1, . . . ,p} are arbitrary indices. 



Considering the cases where the mean vectors \Xt and the covariance matrix 
S take particular values, we will obtain several identities of moments of some 
distributions associated with the Wishart distributions. We shall see that the 
derivations are reduced to enumerating graphs of various types. This is the 
secondary purpose of our paper. 



The Wishart distribution originates with a paper bv lWishart (1928) around 
80 years ago. Since then, it is considered to be a fundamental distribution 
not only in mathematical statistics but al so in other fi elds such as random 
matric es theory and signal processing (e.g., Bai ( 19991) . Maiwald and Kraus 
Despite this, the structure of moments of the Wishart distribution s 
is still an active resea rch to pic. In the central case , iLu and Richards I (|200l[ ). 
Graczvk. et al. I ([20031 ) , and iGraczvk. et al. I ([20051 ) provided formulas for mo- 



ments of the real and complex Wishart distributions. These studies are based 
on expansions of the moment generating functions of the central Wishart dis- 
tributions det(I — 20 ■ The graph presentations of moments have also 
been discussed in these studies. Prior to these studies, it was known that 
terms in the expansio n of det(J — Y)~ a around Y = have some combina- 



torial structures (e.g., IVere- Jones I (|1988h ). which are closely related t o the 
problem of Wishart moment. More recently, iLetac and Massam I ([20081 ) pro- 
vided a method to calculate moments of the noncentral Wishart distribution 
by combining expansions of the moment generating function det(J — 2QE)~? 
and the "noncentral part" e^d^esy^A - m q 



The outline of this paper is as follows. In Section [21 we treat the real non- 
central Wishart matrices. Formulas for the general terms of moments of the 
Wishart distributions are g iven in terms of undirected graphs. This is an ex- 
tension of lTakemura I (|1991l ). who treated the central case. Then, by letting the 
mean vectors /it and the covariance matrix S have particular kinds of struc- 
tures, we ob tain exp l icit fo rmulas for moments of the noncentral chi-square 
distribution, Kibble I (|l94ll) 's bivariate chi-square (gamma) distribution, and 
the 2x2 central Wishart distribution with 2J = I. Noting a formal correspon- 
dence between the moments of the noncentral chi-square distributions and 
the Laguerre polynomials, we will show that the coefficients of the Laguerre 
polynomials have a combinatorial interpretation. 



In Section [31 we treat the case of the noncentral complex Wishart matri- 
ces. Major parts of discussions are parallel to the real case. One remarkable 
difference is that moments in the complex case are not described in terms of 
undirected graphs but directed graphs. 
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2 Moments of the real noncentral Wishart distribution 

2.1 A graph presentation 

In this subsection, we provide a graph presentation formula for moments of 
the real noncentral Wishart distri butions of g e neral degrees. Our results are 



generalizations of Theorem 4.3 of iTakemura I (|1991r) where the central real 



Wishart matrices are treated. Our basic tool is the following formula for mo- 
ments of Gaussian random vectors. T his is just a momen t-cumulant relation 



in the Gaussian case. For the proof, see lMcCullagh I (|1987l) . In the central case 



H = 0, this is sometimes referred to as the Wick formula. 

Lemma 1 (Moment of the real normal distribution) Let X = (xt) be a 
Gaussian random vector with mean /j = (fxi), and covariance matrix E = (<7jy). 
Then, 

E \x\ X2 ' ' ' X 7 i\ — ^ 0~i 1 1 2 0*i2m— 1 M*2m+1 ' ' ' f^in ' 

where the summation is taken over all partitions of n indices {1, 2, ... ,71} into 
unordered m pairs and n — 2m singletons 

[Hi hji • • • 7 (*2m-l, *2m), (*2m+l), ■ • ■ , (in)' 

Remark 1 Although Lemma [T] just states an expression for E\x\Xi ■ ■ ■ x n ], it 
indeed gives general forms of the moments E[x a Xb • ■ ■ x c ] by considering a 
degenerate case. For example, we have E[x\x\] — E\xiX2X^\, where 

0~\2 0~i2 S 
0~21 0~22 0~22 
\&21 O22 0~22, 

Throughout the paper, we will use this degeneracy argument many times. 

Let X t = (xu) (t = 1, . . . , 2/) be independent Gaussian random vectors 
with mean /i t and covariance matrix S. Let W — be a Wishart ma- 

trix made of Xt's as in (JXJ) . In the following, we give a formula for the mo- 
ment E[w a bWcd ■ ■ ■ w e f] with a, b, c, d, . . . , e, / arbitrary indices. By applying 
the degeneracy argument again, we can restrict our attention to the moment 
^[1^12^34 • • • W2n-i,2n] without loss of generality. For example E\w\\w\^\ = 
E[wi2 5)34^56], where ~ W 6 (v, (<7y), (%)), 



Xl\ 






X2 






x 3 J 




\t*2 



'(a n ,5 n ), i,j € {1,2,3,5}, 

(oi2,«i 2 ), i€{l,2,3,5},je{4,6}, 

(o-ai,fei), i e {4,6}, je {1,2,3,5}, 

.(022, 822), «,j'e{4,6}. 

Let V = {1,2,. . . ,2n — 1, 2n} be the set of indices appearing in the ar- 
gument of the expectation E[wi2 ■ ■ ■ W2n-i,2n]- I n the following, we consider 
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an undirected graph whose vertices are the elements of V. First consider an 
undirected graph Go = (V,E ) with the edges 

£J = {(l,2),...,(2n-l,2n)}. 

For each partition of {1, 2, . . . , In— 1, 2n} into m pairs and 2n — 2m singletons, 

(it, i 2 ), ■ • • , («2m-l, «2m), («2m+l), • ■ • , («2n), (5) 

we define a set of undirected edges 

■E = {(h,i2), ■ • ■ , (i2m-li*2m)}- 

By adding the edges of E to Go, we have a graph 

G=(V,E UE). (6) 

Each connected component of G is classified as a "cycle" (a path without 
terminals) and a "chain" (a path with two terminals). For the partition |0, 
the number of chains is n — m. The number of cycles of G is denoted by 
len(G). Note that lcn(G) < m. Let (ji, j 2 ), • • ■ , (j2n-2m-i, 32n-2m) be pairs of 
two terminal vertices of n — m chains of G, and let 

E = \{jl,h)l •••) (j2n-2m-l, J2n-2m)}- 

Using these notations, the general form for the moments is given below. 

Theorem 1 (Moment of the real noncentral Wishart distribution) 

Let (wij) ~ W(u, (<7y), (Sij)). Then, 

E[w l2 ■ ■ ■ w 2n ^ 2n ] = J2 v ien{G) ° E S 6 , (7) 



where 



(i,i')e£ 

TAe summation is taken over all partitions of {1, 2, . . . , 2n} of the form 

Example 1 Consider the evaluation of the moment E[wi 2 w^iW^]. Then, V = 
{1,2,3,4,5,6} and E = {(1, 2), (3, 4), (5, 6)}. There are 76 partitions of V 
into pairs and singletons. Figure Q] is the graph G = (V, Eq U E) for E = 
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1 
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Fig. 1 Graph G = (V,Eo U E) (Eq: solid line, E: dotted line) presenting the term 
v^o i60"25<534 (n = 6, m = 2, len(G) = 1). 



{(1, 6), (2, 5)} (.E = {(3, 4)}). Summing up 76 possibilities, we have the follow- 
ing: 

E[wi 2 w 3 4W 56 ] ^iy 3 <J 12 cr 34 a 5 e + v 2 o 2 zo li a^\ + f<7 2 304 5 <7i 6 [8] 

+ I/ 2 0-12034 (5 56 [3] + ^O-230-i4<W6] + ^120-45(536 [12] 

+ 0-230-45(516 [24] 
+ Wi 2 (334(556[3] + 023(>14<556[12] 
+ 812834655. 

Here, [n] means that there are n terms of similar form. 

Proof For i = 1, ... ,71, let e(i) = [(i + l)/2] (the integer part of (i + l)/2). 
Noting that toy = X^t=i x ti^tj, and from Lemma [TJ we have 

-E[wi2 • ■•1W2n-l,2n] 

= "•• /Z E [ X ti,l%ti>2 • • • X tn .2n-lX tn ,2n\ 

f 1 = l t„=l 

= ^ ' ' '^Z^Ned^ l a teW,2 • • • a; t c (2„-i),2r l -l£t e(2 , l) ,2n] 
ti t„ 

= IZIZ " ' C ° V ( X Un)^ ' X * = (. 2 )^) • • • C0V ( 2; t c ( l2m _ 1 )^ m -l » a! t«(i am ),i2m) 

E ti t n 

X S I a:t e(» 2m+1 )^2 m+ l] ' ' ' ^[ X te(, 2 „),i2„]- (8) 

Since . . . , i„} = V, the indices ii, . . . , i n can be divided into connected 
components of the graph G. Let {ji , . . . , j2fe} be a set of vertices of a connected 
component. Then, as we already pointed out, it forms either a chain 

{jl,h), {32, js), (j2fc-2,j2*-l), C?2fe-1, J2fe), 

or a cycle 

(jl,32), {32,33), • ■ • , (32k-l,32k), C?2fc, jl), 
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and in both cases 

(iiiia), (is, 34), • • • , (j2fc-i)i2fc) e £b. 

Since the running indices ti, . . . ,t n correspond to n edges of i?o> and e(ji) = 
e(j2),e(j 3 ) = e(j 4 ), . . . ,e(j 2 fe-i) = e(j 2 fe), the argument of the summation 
in ((5J) is written as a product of terms of the form 

H " ■ X! E i x tuh] Cav ( x ti,h> x t2,j a ) ■ ■ ■ 

tl tk 

x Covixt^^j^^xt^j^jElxt^^} (9) 

in the chain case, or 

^■■•^Cov(xt lj2 ,x 42j3 )--- 

ti t h 

x Cov(a: tfc _ li:7 - 2ifc _ 2) a;t fc) j 2J ,_ 1 )Cov(a;t fc ,i 2 , > ,a;t 1 ,j 1 ) (10) 
in the cycle case. Here, we used a reindexing 

tl :=i e (ji) = *e(ia)i • • • >*fc := *e(j 2 *-i) = ^e(j 2k )- 

Noting that Cov^^Ey) = l{ s =t}Oy and ELi ^I 1 "]^! 1 *:] = we see 
that © = <r h j 3 • • • aj^j^Sj^ and JTDJ = va hh ■ ■ ■ a hk _^ hh _^a j2kh . 
This completes the proof. □ 



2.2 Enumeration of undirected graphs 

In this subsection, we will enumerate the graphs G = (V, Eq U E) defined in 
(|6]) under the condition that the number / = len(G) of cycles and the number 
m of edges of E are given. Let fi, m ,n be the number of such graphs. 

Consider a degenerate noncentral Wishart matrix W — (wij) such that 
(jjj = 1 and Sij = S. This happens when every component of X t making up W 
takes the same value with probability one. Accordingly, all elements of W take 
the same value w, say, with probability one. In this setting, ([7]) in Theorem 
[T] is reduced to a moment formula for the distribution of w, the noncentral 
chi-square distribution xt(^) with v degrees of freedom and the noncentrality 
parameter S. Using the coefficient fi >m ,n, the nth moment of w is given as 
follows. 

n 

m=0 Z>0 

The coefficient fi iT n,n satisfies the following recurrence formula. 
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Lemma 2 

fl,m,n = 2(2n — m — l)/i,m-l,n-l + //-l.m-l.n-l + /i,m,n-lj (12) 
im'i/j boundary conditions 

f Ji a = o) 

(i o „ = < for n > 1, 1<5 

J ' ' [0 (J>1) " 



f -J° (/ = 0) ' run 

a = 1) _ (14) 

Proof In the following, we sometimes refer to an edge from E$ as a "solid line" 
edge and an edge from E as a "dashed line" edge as shown in Figure [TJ 

The connected components of G are classified as cycles and chains. In each 
cycle, the number of solid line edges is equal to the number of dashed line 
edges. In each chain, the number of solid line edges is one more than the 
number of dashed line edges. Thus, the number of connected chains is n — m, 
the difference between the number of solid line edges and the number of dashed 
line edges. 

Consider a graph G made by removing two vertices 2n — 1 and 2n, and all 
(solid line and dashed line) edges connecting to these two vertices. Note that 
the (solid line) edge (2n — 1, 2n) is deleted. 

One of the following will meet: The edge (2n — 1, 2n) is contained in (i) a 
cycle with 4 or more edges, or a chain with 3 or more edges (the edges (a, b), 
(c, d), or (fc, I) in Figure[2]); (ii) a cycle with 2 edges ((e, /) in Figure[2]); 

(iii) a chain consisting of a edge {{g 7 h) in Figure [2]). 

Case (i). In the graph G' made by removing two vertices In — 1 and 2n, 
and all connected edges, there are n — 1 solid line edges. Since one dashed line 
edge is removed together, there are m — 1 dashed line edges. The number of 
chains still remains (n — 1) — (m — 1) = n — m. To the graph G", consider 
adding the edge (2n — 1, 2n) again. There are n — 1 + (n — m) = 2n — m — 1 
places where (2n — 1, 2n) can be inserted, and considering the direction of the 
inserted edge, there are 2(2n — m — 1) ways in which the insertion can be 
done. By this operation, the number of dashed lines increases by 1, whereas 
the number of cycles is invariant. The contribution of the number of graphs 
made by this operation to fi <m ,n is 

2(2n -m - l)/i >T n-l,n-l ■ 

Case (ii). Consider adding the edge (2n — l,2n) to the graph G' again 
to make a cycle with the edge (2n — l,2n) and a dashed line edge. By this 
operation, both the number of dashed lines and the number of cycles increases 
by 1. The contribution of the number of graphs made by this operation to 

fl,m,n IS 

fl— l,m— l,ro— !• 
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a 

b I 



a^9 



Fig. 2 A figure for the proof of Lemma [2] 



f 9 " 



Case (iii). Consider adding the edge (2n — 1, 2n) to the graph G again to 
make a chain consisting of one edge (2n — 1, 2n). By this operation, both the 
number of dashed lines and the number of cycles are invariant. The contribu- 
tion of the number of graphs made by this operation to fi. m ,n is 



fi 



l,m,n—l ■ 



Summing up the three cases (i), (ii), and (iii), we obtain the recurrence 
formula d. □ 



Theorem 2 The generating function of //,,„,„ with respect to the number I of 
cycles is given by 



l>0 

n 
m 



$m,,nty) = y"V/i,m,n 
m 

TT(i/ + 2(n-i)) (0<m<n, n > 1). (15) 
«=i 

Here, we use a convention Yii=i = 1- 

Proof Noting that /_i, m , n = 0, the generating function $ m , n (i>) = J2i>o ul fl,m,r. 
has to satisfy 

<£ m ,„ = 2(2n m l)<P m _ 1>n _i + V$ m ^ x ,n-\ + #m,n-l- (16) 

To solve the recurrence formula (|16p. consider the boundary conditions. 
From ([TU|) . we have 

M") = l («>!)• (17) 
Moreover, from (fT^j) . we have 

$1,1(1/) = ZA 
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Furthermore, since <? ni „_i = 0, 

= 2(n - l)#„_ 1)tl _! + ^„_ 1>n _ a = (1/ + 2n - 2)<5 n _ 1 ,„_i = • • • 

= + 2(n - i)) (n>l). (18) 

i=l 

The recurrence formula (|16p combined with the boundary conditions (|T7| and 
([T8)l determines $ m , n for all m and n. 

In the following, we see that $ mi „(i/) in (|15l) is actually the solution for the 
recurrence formula (fTfj|) . The boundary conditions (|17p and (fT5|) are satisfied. 
We only have to make sure that (fT5)l really satisfies (TTB)) . Indeed, 

" ^m,n-l = ( ) JJ^ + 2n ~ 2i ) ~ ( ) IE" + ^-2-2l) 



m 

n-l\ 1 



i=l v ' i=l 

m— 1 



m — 1J m 



rr(f + 2n-2-2i) 

i=l 

x + 2n - 2) - (n - m)(V + 2n - 2 - 2m)} 
2(2n - m - l)$ m _i >n _i + 2/$ m _i >n _i. 



□ 



Corollary 1 

^ron(l) = | )(2n- l)(2ra-3)---(2n-2m + l) 

\TO J 

2n \(2m-l)\\ 
is the number of undirected graphs G, and 

*m,n(0) = f n J (2n - 2)(2n - 4) ••• (2n - 2m) 



2 m n!(n- 1)! 



to! (n — to)!(ti — to — 1)! 
zs the number of undirected graphs G without cycles. 

Remark 2 Nonnegative integers s n (m,l) defined by a generating function 

tci m 

£Va n (m,Z) =Y[(v + n-i) (19) 

1=0 i=l 

are called the noncentral Stirling numbers of the first kind ( Koutras (1982)). 
Since 

5>7*,m,n= ( n )2 m f[(u/2 + n-i)= Q2 m 5>/2)'.s„(m,0, 

u>0 i=l l>0 
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wc have 



\m — l 



s n (m,l). 



2.3 Moments of the noncentral chi-square distribution and the Laguerre 
polynomial 

As stated in the beginning of the previous subsection, the moment of the 
noncentral chi-square distribution is described with the coefficient fi, m ,n- 

In view of (jlll) and Theorem [2j the nth moment of w ~ xt{$)i the noncen- 
tral chi-square distribution with v degrees of freedom and the noncentrality 
parameter S, is written as 



E\w n 



m=0 l>0 
n 

n / \ m 

= E(")n("+ 2 (»-o^ 



m=0 



(20) 



This is a well- known expression for t he moment of noncentral chi-square dis- 
tribution (e.g., Johnson, et al. ( 19951 ) V 

Remark 3 iKoutras" (1982) pointed out that moments of some noncentral dis- 
tributions are described with the noncentral Stirling numbers of the first kind. 

The moment generating function of the noncentral chi-square distribution 



This can be obtained by letting = t, £ =1, A = S with ©. 

The Laguerre polynomials, the orthogonal polynomial systems on (0, oo) 
with respect to the gamma weight functions, are defined as 

L^(x) = 2 n — f^(x;a) a J&\^ 1) (* > 0), 



where 



f^{x;a)=x v ' 2 - x a- v ' 2 e- x '^ 



(e.g.. lMorris I ([19821) ). From this definition, we immediately get the generating 
function of the Laguerre polynomial as 



X)(-i) n - T 4" ) (^) = (i-2t)- 



"/2g-0/2)((l-2t)- 



-1) 



n=0 
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which has formal coincidence with the moment generating function of the chi- 
square distribution with v degrees of freedom and the noncentrality parameter 
— x. Therefore, we have an expression for the Laguerre polynomials 



l%\x) = (-1)- ^^^^(-^r 

n / \ m 



m=0 



This gives a combinatorial interpretation for the coefficients of the Laguerre 
polynomials. 

This type of combinatorial interpretation for Hermite polynomials is widely 
known. Applying a degenerate multivariate distribution AT((/i,), with 
fii = /i, Oij = a 2 to Lemma [TJ we see that the nth moment of the normal 
distribution X ~ N(fj,, a 2 ) is 

[n/2] 

prynl \ 2m n— 2m 

m—0 

where 

a mn =[ (2m - 1) = - r 

' V 2m / (n-2m)!2™m! 

is the number of partitions of an n-member set into (unordered) m pairs and 
n — 2m singletons. The Hermite polynomials are defined by 

H n (x) = { -i r £-e-* 2 /*/e-^. 
The generating function of H n (x) is 

00 j-n 
n=0 n ' 



which coincides with the moment generating function e Alt+cr * I 2 of the nor- 
mal distribution N(fx,a ) with /x and cr 2 replaced by a; and —1, respectively. 
Therefore, the nth Hermit e polynomial is the n moment of the norma l distr i- 
bution N(x,-l) formally (iMcCullagh I (|l987l) . iKuriki and Takemura I (|l996h . 



Withers and Nadaraiahl ( 20061 )1 . and hence 



[n/2] 
m=0 
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2.4 Moments of bivariate chi-square distribution 

In this subsection, we give an explicit expression for the moment of bivariate 
chi-square distributions as the second application of the graph presentations 
for the moments of the noncentral Wishart distribution. 

There are several proposals for defining bivariate chi-square (gam ma) dis- 
tribut ions. The distribution we will discuss here is that given by iKibble 



(|194l[ ). Kibble's bivariate chi-square distribution is defined as the distribu- 
tion of the diagonal elements (wn,w 22 ) of a 2 x 2 central Wishart distribution 
(wij) ~ W 2 {v,S) with 

Substituting — diag(*n, £22) and A = into @, we have the moment 
generating function 

p\ f ,t 11 w 11 +t 2 2w 22] _ j . A - 2tu -2/otii \ / 

= (1 - 2t n - 2*22 + 4*11*22 - 4p 2 t n t 22 )- V/2 . (21) 

The next theorem gives an expression for the moments of general degrees 
E\w\iW2 2 \- We evaluate this as the moment 

E[w\2 ■ ■ ■ W 2 b-1, 26*026+1, 26+2^2(6+c)-l,2(6+c)]) 

where ~ W 2 ( b+C )(v, (o^)) with 

jl <2b or i,j > 26 + 1), 

I p (otherwise) . 

Theorem 3 (Moment of the bivariate chi-square distribution) Let b 

and c be nonnegative integers. Then, 

min(&,c) 

E[ W w 22] = £ ^ ft. a)1 ( -: 0)1 d n(»+g(«-o) 

a=0 v ' v ; z=l 



6— a ■ " 

X 

i=l 



n^+2(6- i ))n^+ 2 ( c - i ))- 



Remark 4 iNadaraiah and Kotzl (|20Q6l) derived an expression for Elw^w^} in- 
cluding the Jacobi polynomials. Their derivation is a use of a series of identities 
of special functions, which are totally different from our combinatorial proof 
given below. 

Proof Assume that b < c without loss of generality. Let Go = (Vi UVi, EiUE 2 ) 
be a union of two undirected graphs (Vi, Ei) (i — 1, 2), where 

VI = {1, 2, ... ,26 - 1, 26}, Eh = {(1, 2), . . . , (26 - 1, 26)} 
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and 



V 2 = {2b + 1, 2b + 2, . . . , 2{b + c) - 1, 2(6 + c)}, 
E 2 = {(2b + 1, 26 + 2), . . . , (2(6 + c) - 1, 2(6 + c))}. 



By forming 6 + c pairs from 2(6 + c) vertices of V = V± U V%, we add 6 + c 
edges joining two vertices of each pair to the graph Go to make G. Let ft./, a ,b, c 
be the number of resulting graphs G such that the number of cycles is I and 
the number of edges joining an element of V\ and an element of Vi (the number 
of pairs consisting of an element of V\ and an element of V-i) is 2a. Then, the 
moment that we want to evaluate is represented as 



We divide the process of adding b + c edges into three steps below (Figure [3]) . 

Step (i). Choose 2c — 2a from the 2c vertices of V2, and form c — a pairs 
from the 2c — 2a vertices. Add c — a edges joining two vertices of each pair to 
the graph Gq . The number of resulting graphs having I' cycles is 



Note that 2a vertices not chosen are the terminal vertices of a chains. 

Step (ii). Form 6 pairs from the 26 vertices of V\, and add 6 edges joining 
two vertices of each pair. The number of resulting graphs having I" cycles is 



Step (hi). Choose a edges from the 6 edges added in step (ii), and make 
"cuts" at the middle of each edge, and have the a chains generated in step (i) 
fit in at the a cut points. Note that this operation does not alter the number 
of cycles. Since the a chains have directions, the number of ways to perform 
this operation is 



Summing up (i), (ii), and (hi), we get 



^v l p 2a h UaAc . 



I . a 



fl',c 



fl",b,b- 




Therefore. 



^ v l p 2a hi, a ,b,c 



Hp- 



2 a 6! 




/ . a 



a 



(&-<.)! 



s^(^C-„)n<'' +2 ( c - ! »n('' +2 < t - i »»- 



□ 
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I I 
I I 

I * 

I f 



Fig. 3 Figures for the proof of Theorem [3] (Left: Steps (i), (ii), Right: Step (iii); O : Vertex 
of Vi, • : Vertex of V 2 ). 



2.5 Moments of a 2 x 2 real Wishart distribution 



As the third example, we give an explicit expression for the moments of a 2 x 2 
central Wishart distribution with the parameter S = I. Let (t£>y) ~ W^v,!). 

Substituting = ( ,^ 11 /r , ^ | and A = into ©, we have the moment 

\I12/2 t22 / 

generating function 

- Of -f x ~" /2 



^|gtl2«'12+tllU>ll+t22tU22 j 



-tl2 1 — 2£ 22/ 

(1 - 2i u - 2t 22 + 4< n < 22 - i? 2 )~ l//2 . (22) 



We first show that 

E\w\ 2 w\ x w\^ =0 for a odd. 



Let X t be Gaussian random vectors making up the Wishart matrix. Since 
E[X t ] = 0, the distribution of X t is invariant under the change of the sign of 
the first coordinate. On the other hand, this change causes (wi2, Wu, 11*22) ^ 
(— w%2, W11, 11)22)- This implies that 

E[w a l2 w\ lW c 22 ] = E[(~w 12 ) a w\ lW c 22 ] 
= (-l) a E[w a 12 w b llW c 22 }. 



Unless a is even, the left- and right-hand sides become 0. 

In the following, we will derive the moment £/[w?2 w xi w 22\> w here a, b, and 
c are nonnegative integers. Although some methods to calculate this quantity 
are already known (e.g., Remark [5]), we demonstrate that our combinatoric 
approach does get the same results. 

Let 



= {1, 3, . . . , 4a - 1}, = {4a + 1, 4a + 2, . . . , 4a + 26}, 
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and 

V 2 = {2, 4,..., 4a}, V 2 + = {4a+ 26+ 1,4a + 26 + 2,..., 4a + 26 + 2c}. 
Then, 

^[^12^11^22] = E[w 12 W 34: ■ ■ ■ «'4o+26+2c-l,4o+26+2c], 

where (wiy) — W± a+2 b+2c{v-, (Vy)) with 

Ji e V+, or i, j e V+, or i e e V 2 °, or i e V 2 °,j € 

(7V j = < 

10 (otherwise). 

Theorem 4 (Moment of the 2x2 real Wishart distribution) Let a, b, 

and c be nonnegative integers. Then, 

a b c 

E[w\^w\ lW c 22 ] = {2a-l)l\]J(u+2(a-i))Y[(u+2(a+b-i))Y[(u+2(a+c-i)). 

i—l i—1 i—1 

Proof Let V\ — V® U and V2 = V 2 ° U V 2 ' . First, define an undirected graph 
G = (V, E Q ) with 4a + 26 + 2c vertices V = Vi U V 2 and 2a + 6 + c edges 

£ = {(1,2), (3, 4),..., (4a + 26 + 2c- 1,4a + 26 + 2c)}. 

Then, consider the addition of additional 2a + 6 + c edges to the graph Go to 
make G such that no edges joining V\ and V2 are added. Let h^ a ,b,c be the 
number of resulting graphs G such that G has I cycles. Then, 

l>0 

We divide the process of adding 2a + 6 + c edges to Go into three steps (i), 
(ii), and (iii) below (Figure @|. 

Step (i). Let < a' < min(a, c). Choose 2c — 2a' elements from the 2c 
vertices of V 2 , and form c — a! pairs from them. Add c — a! edges defined by 
the c — o! pairs to the graph Go- According to this operation, a! chains are 
newly generated. The number of graphs with V cycles is 

fl',c—a',e' 

Step (ii). Form a pairs from the 2a vertices of V 2 . The number of ways in 
which this pairing can be done is 

(2a-l)!! = M. 
v ; 2 Q a! 

Choose a' pairs from the a pairs, and assign each pair to each of the a' chains 
generated in step (i). Connect a vertex of the pair to one terminal vertex of 
the chain using a new edge, and connect the other vertex of the pair to the 
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other terminal vertex of the chain using another (new) edge. (Add 2a' edges 
in total.) The number of the correspondences is 



2 a a' 

2a(2a - 2) ■ ■ ■ (2a - 2o') 



(a-a'y: 

For the remaining a — a' pairs, connect two vertices of each pair using a new 
edge. (Add a — a' edges in total.) 

The number of edges added in steps (i) and (ii) is (c — a') + 2a' + (a — a') = 
a + c. The number of cycles is I'. 

In steps (i) and (ii), summing the number of ways for < a' < min(a,c) 
yields 

min(a,c) ^ 



:,a= J2 /l',c-a',cX(2a-l)!l 



X 



(a -a')! 

a'=0 v ; 

The coefficient ei>, c ,a can be combinatorially interpreted as follows: Let G 2 = 
{V 2 ,E 2 ) with E 2 = '{(4a+26+l,4o + 26+2),...,(4a+26+2c-l,4a+26+2c)} 
be an undirected graph, and add a + c edges by forming a + c pairs from the 
2a + 2c vertices V 2 . Then, a chains are newly generated, ey jCja is the number 
of resulting graphs having V cycles, and the terminal vertices of the a chains 
are elements of V 2 °- 

Step (iii). For all 2a + 2b vertices of V\, form a + b pairs and connect the 
vertices of each pair with a new edge. (Add a + b edges in total.) According 
to step (ii), the 2a vertices of V\ have already been divided into a pairs, and 
the vertices of each pair have been connected with an edge. In step (iii), the 
number of graphs adding new I" cycles is fi» , a +b,a+b- 

Summarizing (i), (ii), and (iii), we get 



hl. a ,b,c — £l' ,c,afl" ,a+b,a+b, 

l'+l"=l 

and hence, 

^^V l hi^ a ^b,c = yl e l'^M ^ V 1 fl" t a+b,a+b- (23) 
i>0 l'>0 l">0 

For a nonnegative integer n, write 

'a\ a(a — 1) • • • (a — n + 1) 



n 



7! 
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Fig. 4 Figures for the proof of Theorem [4] (Left: Steps (i), (ii), Right: Step (iii); □ : Vertex 
of Vj , O : Vertex of V+, ■ : Vertex of V 2 °, • : Vertex of V+). 



Then, the generating function of the coefficient ey iCj(1 with respect to the num- 
ber of cycles V is 

5>V C;0 = (2a-l)!! £ E^/^-V7^M 

i'>0 a'=Q l'>0 ^ '* 



min(a,c) 



c— a 



<*-»» e t-jn^*-^ 

a'=0 V 7 i=l V ; 



min(a,c) 



)» e (,_ e j^r*r i v«-^"' 



a'=0 

min(a.e) 



c—a'J \ c — a' J (a — a')! 



(2a-l)!!2 c c! ^ 



a'=0 



a\ fv/2 + c- 1 
a' J I c - a' 



= (2 a ~l)!!2^!(^ 2 + ; +C - 1 

C 

= (2o- 1)!! JJ(i/ + 2(o + c-i)). 
i=i 

The fifth equality above is known as a convolution identity for two binomial 
coefficients. Substituting this into (|2"5)l completes the proof. □ 



Remark 5 Theorem 3] can also be shown by the following geometric consider- 
ation. Let X* — ( 41 ) be Gaussian random vectors making up the 2x2 

\Xt2j 

Wishart matrix (wtj) ~ W2[v,I)- Write X^ = (xu, . . . , x V i)' (i = 1,2), 
where ' denotes the matrix transposition. Since the ^-dimensional distribu- 
tion of X(f) is invariant under the orthogonal transformation preserving the 

norm ||-X7j)|| = ,/X',.sXu\, four quantities ||-X7j)||, _X7j)/||.X"(j)|| for i = 1,2 are 
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independently distributed. Thus, 

E[wl a 2 w b n w c 22 ] = J B[(X ( ' 1) X (2) ) 2a ||X (1) || 2fc ||X (2) || 2c ] 

= i?[r 2a ]S[||X (1) || 2 ( a +'')]i?[||X (2) || 2 ( a + c )], 

where r — X^X(2)/(||X( 1 )||||X( 2 )||). This calculation can be completed by 
noting that HX^H 2 ~ xt an d 7,2 ~ ^-), the beta distribution. 



3 Moments of the noncentral complex Wishart distribution 

3.1 Preliminaries on the complex normal distribution 

In this section, we will deal with the complex noncentral Wishart matrices. 
Major parts of the discussion are parallel to the real case. One remarkable 
difference is that the moments in the complex case are described in terms 
of directed graphs, whereas those in the real cases are described in terms of 
undirected graphs. 

We begin by summarizing some preliminaries on the complex normal dis- 
tribution and the complex Wishart distribution. Let Z = (zl) be a complex 
conjugate of Z = (Zi)- The following lemma is a complex version of Lemma [T] 

Lemma 3 (Moment of the complex normal distribution) Let Z = 

(z,^ <~ CN(ji, S), and let i = n + i, i = 1, . . . , n' . Then, 



E[z x --- z n z T ... ZT7 ]=^2 a iih ■ ■ ■ a imj, 



where the summation is over possible pairing ji), (i2>j2)> • • • (im,jm)} such 
that 

{h,..., i m } C {1,2,..., n} and {j x , . . . ,j m } C {T, . . . , ri} 
(i.e., matching). The other indices are 

{i m+ i,...,i n } = {l,...,n}\{ii,...,j fc }, 

{jm+l, ■ ■ -Jn'} = {T, ...,»'} \ {jl, ■ . -,jk}- 

Proof Write Z = X + y/^lY, fi = Z + y/^ln, £ = 2(A + ^/-LB). Let 9 = (^), 
if = (ipi) be parameter column vectors. Because of 

X\ . N ((t\ (A-B 



and 



Y J \\vj ' \ B A 
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the moment generating function of (Z, Z) is obtained as 
E[exp{9'Z + V 'Z}] = cxp{ (« v ') J + ± (9' /) J (g "/) (*) } 

-4sV)(£)4(»V)(^)C)} 

= expjflV + yp'/J -(- ~(0' Sip + tp'£9)} 
= exp{6*'/i + ip~p + 9' Uip}. 

From this, we have the joint cumulants of (Z, Z) as 



/'l 


(n 


= 1, 


n' 


-0), 


Mr 


(n 


= o, 


n' 


= 1), 


°1T 


(n 


= 1, 


n' 


= 1), 



[0 (otherwise). 

Lemma |3] is the moment-cumulant relation for this particular cumulants. □ 



3.2 A graph presentation 

Let Z t — (z t i) (t = 1, . . . , v) be independent complex Gaussian random vectors 
with mean fj, t and covariance matrix S. A complex Wishart matrix W = (wij) 
is constructed from Z t as given in In this subsection, we give a formula 
for the moment E[w a bWcd ■ ■ ■ w e f] with a, b, c, d, . . . , e, / arbitrary indices. By 
considering the degenerate case again, without loss of generality, we only have 
to treat the moment E[w 1 jiv 2 2 1 ■ ■ w-an] with i = n + i, i = 1, . . . , n. 

Let V = {1,2, ... ,n} be the set of indices appearing in the expectation 
that we want to evaluate. In the following, we consider a directed graph whose 
vertices are the elements of V. Choose a subset V\ of V such that its cardinality 
is \Vi\ — m, and consider an injection tt : V\ — > V. The map it defines a set of 
directed pairs 

E = {(i,n(i))\ieV 1 }. 

We regard the pair (V, E) as a directed graph G, where V is the set of vertices, 
and E is the set of directed edges. Note that E and the pair (Vi, tt) have one- 
to-one correspondence. 

As in the undirected case, for a given G, every connected component is clas- 
sified as a "cycle" (a directed path without terminals) and a "chain" (a directed 
path with a starting terminal and an ending terminal). For the map tt, the num- 
ber of chains is n—m, where n = \V\, m = \V\\. The number of cycles of G is de- 
noted by len(G). Note that len(G) < m. Let (ji,kx), . . . , (j n -m, k n -m) € VxV 
be directed pairs of ending and starting terminal vertices of n — m chains of 
G, and let 

E = • ■ • , (jn-m) kn-m)}- 
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o 



3 



Fig. 5 Directed graph G = (V, E) presenting the term f 1 o" 1 xf 23^32 = 3, m = 2, len(G) 
1). 



Using the notations above, we give the general form for the moments as 
follows. 



Theorem 5 (Moment of the complex noncentral Wishart distribu- 
tion) Let (wij) ~ CW{v, ((Jij), (5ij)), and let i = i + n, i = 1, . . . , n. Then, 

E[w lT ■ ■ ■ Wrm ] = J2 v l ™ [G) v E ^, (24) 



where 



IT u ii> — u ii7r(ii) " ' "i m ir(i m )> 



» ,: n h 



33 3\3\ 3n- m 3„ 



The summation is taken over all possibilities of V\ = {i±,...,i m } C 

{1, . . . , n}, and injections n : Vi — >• V. 

Example 2 Consider the evaluation of the moment E[w 1 jw 2 2Wgg\. Then, V = 
{1,2,3}. There are 34 injections from subsets Vi C V to V. Figure [5] is the 
graph G = (V,E) for E = {(1, 1), (2, 3)} (E = {(3,2)}). Summing up 34 
possibilities, we have the following: 

£KT W 22 W 33] = i/3cr lT Cr 22 Cr 33 + ^ °12 I 3 ] + "°'l5 CT 23 £7 3T [ 2 ] 

+ o-^o-^S^] + ^CT 1 2CT 2T (5 3 3[3] + W lT (T 2 3(532[6] + (T^a^S^ [6] 

+ (5 lT (5 22%- 



Here, [n] means that there are n terms of similar form. 
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Proof Note that Wij — J^t=i z ti z tj- In view of Lemma O we have 

E[w lT ■ ■ ■ W n n] 

v v 

= / J "•• / ] ^[ z ti,lZ tl j- ■ ■ Zt n>n Zt n ,n] 
t 1= l t„=l 

= S * * * S Cum ( z *n A ' z t„ (il) ,¥(iT)) ' ' ' Cum K m A» . z t, (lm) ,^)) 

B ti t„ 

x ' ' ' E \ z uM E K, F=] • ••£[%7lH. (25) 

*m+l m + 1 " 

where V\={i\,.. .,i m }, 

{i m+1 ,...,i n } =V\Vi, {C+u--->0 = V"\tt(7i). 

Since . . . , i„} = V, the indices i\,...,i n can be divided into connected 
components of the graph G. A connected component having vertices j%, . . . , 
forms either a directed chain 

(h,h), (32,33), (jfe-2, jfc-i), (jfc-i, jfc) 
or a directed cycle 

(iijj'2), (32,33), (jk-l,3k), (jk,h), 

where 7r(j'j) = ji+i (and 7r(j'fc) = ji in the cycle case). Since the running n 
indices ti, . . . , t n correspond to n vertices of V, the argument of the summation 
J2e 1X1 i s a product of terms of the form 

^[^~jr] Cum ( z * 31 a . %J;) Cmn ( z h 2 , ) • • • 

x Cum(z t . fc _ i , ; ,- fc _ 1 ,^— ^ )E[z tjk ljk ] (26) 
in the chain case, or 

x Cum^^^^^Cum^j,,^) (27) 

in the cycle case. Noting that 

Cum(z SI ,z^) = Cov(z si , z t j) = l{ s = t }<7- 

and Et=i ^fei]^] = %, we see that @ = ^ - • • • a Jk _ -6^- and 
<G3> - w ll3 jff A5 • • ■ <t, -<7, 5-. This completes the proof. □ 
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3.3 Enumeration of directed graphs 

In this subsection, we evaluate the number of directed graphs appearing in the 
expression (|24[) for the moments of the complex Wishart distribution. 

Let V = {1, . . ., n}. Let V\ be a subset of V such that \V\\ = m. Let ir be 
an injection V\ — > V. As explained in the previous subsection, we can define 
a directed graph G = (V,E) with E = {(i,ir(i)) \ i e Vt}. The connected 
components of G arc cither a directed cycle or a directed chain (the length 
may be 0). Note that the number of chains is n — to. Let gi >m ,n be the number 
of such graphs having I cycles. The coefficient gi, m ,n satisfies the following 
recurrence formula. 

Lemma 4 

gi,m,n — gi-\,m-\,n-l + 9l,m,n-l + (2™ ~ m ~ l)<ft,m-l,n-l (28) 

wi£/i boundary conditions 

9/ o n = s , for n> 1. (29 

y ' ' 10 (/>!) ' " 



jo (/ = o), 

ftA1= \l (/ = !). (30) 

Proof We consider removing the vertex n and the adjacent edges from the 
graph G. There are three types of status about adjacent edges. 

Case (i). 7t(n) = n. In this case, the vertex n is contained in a cycle with 
length 1. Removing the vertex n and the edge (n, ir(n)) yields a graph whose 
values oi l,m and n are one less than those of G. This corresponds to the first 
term in the right-hand side of (|28p. 

Case (ii). Neither 7r(n) nor 7r~ 1 (n) exists. In this case, the vertex n is an 
isolated vertex. Removing the vertex n yields a graph whose I and m are 
invariant and n are one less than that of G. This corresponds to the second 
term in the right-hand side of (j!?8l . 

Case (iii). Otherwise. In this case, the vertex n is contained in a cycle with 
length more than or equal to 2, or contained in a chain with length more than 
or equal to 1. Remove the vertex n and one edge adjacent to the vertex n. 
This manipulation yields a graph whose / is invariant, and whose to and n are 
one less than those of G. Conversely, when we rebuild the graph G from this 
smaller graph, there are (n— l) + {(n— 1) — (to — 1)} = 2n — m— 1 places where 
the vertex n and one edge can be inserted. This corresponds to the third term 
in the right-hand side of l|28p. □ 



Theorem 6 The generating function of the coefficient gi iT n,n with respect to 
the cycle number I, 

l>0 
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is given by 

& m ,n{ v ) = ( " ) Yl^ + n-i) (0<m<n, n > 1). (31) 



Here, we use a convention = 1. 

Proof Because g-i tm ,n = 0, we see 

#m.,n = ^#ro-l,n-l + #m,n-l + ( 2rl - m - l)^ m _i, n _i. (32) 

The boundary conditions $o,n = 1 (ra > 1) due to (|29| and $14 = v due to 
(|3"0]l satisfy (l3Tj) . In addition, since \P n , n -i = 0, 

n 

*n,T. = + n - l)^ n -i,n-i = • ■ ■ = JJ(i/ + n - i) (n > 1). 

i=l 

Hence, it is sufficient to ensure that pip satisfies (|3"2"|) . Indeed, we have 

- «pm,n-i =( ™ ) ri(^+"-*) - ( n m 1 ) fl^+ri- 1 -*) 



■;=i 

m— 1 



71- 1A 1 

m 



1 - n^ +n - i - i ) 

1J m ±A - 

x {n(v + 71 — 1) — [n — m)(y + n — 1 — m)} 
(2n - m - l)^ m -i,„-i + ^ m -i,„_i. 



Remark 6 Comparing (|3"Tj) with (|19jl in Remark [21 we have 



□ 



9l,m,n = }Sn(m,l), 

where s n (m, I) is the noncentral Stirling numbers of the first kind. In particu- 
lar, gi } n,n = s n (n, I) is the Stirling numbers of the first kind. Let V = {1, . . . , 71} 
and let 7r : V — > V (bijection). It is well-known that the Stirling number of 
the first kind s n (n, I) is the number of directed graphs (V, E), E = {(i, n(i)) \ 
i G V"} having I cycles ([Stanley I (|2000h ). 



Corollary 2 



2 

71 \ , v / 71 



= n(n- l)---(n-m + l) = m! 
is i/ie number of directed graphs G, and 

n\(n — 1)! 



«Vn(0)= (n-l)(n-2)...(n-m) 



ml (n — m)!(n — m — 1)! 
is i/ie number of directed graphs G without cycles. 
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3.4 Degenerate cases 

3.4-1 The noncentral chi-square distribution 

As in the real case, we can obtain several identities for moments by assuming 
that the parameters S and A have particular kinds of structures. First, we 
consider the case where S = (cy), cry = 2, and A — (<%), <% = S. Then, 
every element of W has the same value w, say, with probability one, and 
the distribution of w is the noncentral chi-square distribution xiL(^) with 2i/ 
degrees of freedom and the noncentrality parameter S. The nth moment of 

w ~ xL(£) is 

E[w n ] = E[WjJ- ■ ■ W n n} 

n n 

= E Y.^9l,m,n2 m S n - m = *m,n(v)2 m t n ~ m 
m=0 l>0 m=0 
rt y -v m n , \ m 

= E U n(" + » - *) 2W * n ~ ro = E u IK 2 - + 2 (« - w n ~ m 

m=0 ^ ' i=l m=0 ^ ' i=l 

This coincides with the formula ([2"U]) obtained from the real Wishart distribu- 
tion. 



3.4-2 Moments of a bivariate chi-square distribution associated with the 
complex Wishart distribution 

Let (wij) ~ GWiiy, S) be a 2x 2 central complex Wishart matrix. We consider 
a particular structure of the parameter: 




The diagonal elements (ton, W22) are distributed according to a sort of bivari- 
ate chi-square distribution, since the marginal distribution of wn and W22 are 
the chi-square distribution x%v anc ^ they are correlated. At a glance, (u>n, W22) 
has a different distribution from Kibble's distribution since it has a different 
origin. However, from (jj), the moment generating function is 

E[e t^ 11+t22W22] = det A - 2t u -mn \ 

\ -2/9*22 1 - 2*22/ 

= (1 - 2t u - 2* 22 + 4*11*22 - 4^*11*22)^, 

which is equal to the moment generating function (|21j) of Kibble's bivariate 
chi-square distribution with v and p replaced by 2v and respectively. 
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3.4-3 Moments of a 2 x 2 complex Wishart distribution 

Consider a 2 x 2 complex Wishart matrix (wjy) ~ CW 2 {v, I). We first show 
that 

Let Z t be complex Gaussian random variables from which the Wishart matrix 
W is constructed. Since E[Z t ] = 0, its distribution is invariant when the first 
element of Z t is multiplied by y— 1. On the other hand, this manipulation 
causes 

w\i wi 2 \ ( w\\ \f-iwyi 

W21W22J \— \f— lw-2.1 W22 

Therefore, 

E[w a l2 w^w b lW c 22 ] = E[(V=lw 12 ) a (-V=iw 21 r'w b u w c 22 ] 



-a — a 



= V-l" " E[w1 2 W2iW b n w^ 2 

The left- and right-hand sides become unless a = a' . 

Substituting = (Uj)' and A = into fllj), we obtain the moment gener- 
ating function 



— // 



_ 1 _ t 2 2y 

= (1 - <n - t 22 + t u f 2 a - t 12 t 21 )-». (33) 

Let 



6„(it, u, w) = (1 — u — V + uv — w 



) ^\u,v,w) = (J^J <t>v{u,v,w). 



Then, 

m = Mhi,tzt,ti2t2i) = 4 fc) (*u.*aa,0). 

fc>0 

On the other hand, the moment generating function of the 2x2 real Wishart 
distribution is rewritten as 

t 2k 

m = </v/ 2 (2in,2Wi 2 ) = E -fC/U^U, 2*22,0). 

k>0 

Comparing the two functions, we can immediately obtain the moments of the 
2x2 complex Wishart distribution from the results for the real case in Theorem 

m 

Theorem 7 (Moment of the 2x2 complex Wishart distribution) Let 

a, b, and c be nonnegative integers. Then, 

a b c 

E{w b 11 (wi2W2i) a w 2 ] 2 ] = a! Y[(i> + a-i) \[{v + a + b-i) Y[(i> + a + c-i). 
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